EUROPHYSICS LETTERS 
Europhys. Lett., (), pp. () 



Superconducting correlations in ultra-small metallic grains 

Guang-Shan Tian 1 ' 2 , Lei-Han Tang 1 and Qing-Hu Chen 1 ' 3 

1 Department of Physics, Hong Kong Baptist University, Kowloon Tong, Hong Kong 

2 Department of Physics, Peking University, Beijing 100871, China 

3 Department of Physics, Zhejiang University, Hangzhou 310027, China 

(received ; accepted ) 

PACS. 74.20Fg - BCS theory and its development. 

PACS. 71.10— w- Theories and models of many electron systems. 

PACS. 71.24+q - Electronic structure of clusters and nanoparticles. 



Abstract. - To describe the crossover from the bulk BCS superconductivity to a fluctuation- 
dominated regime in ultrasmall metallic grains, new order parameters and correlation functions, 
such as "parity gap" and "pair-mixing correlation function", have been recently introduced. 
In this paper, we discuss the small-grain behaviour of the Penrose-Onsager-Yang off-diagonal 
long-range order (ODLRO) parameter in a pseudo-spin representation. Relations between the 
ODLRO parameter and those mentioned above are established through analytical and numerical 
calculations. 



The hallmarks of superconductivity — zero resistance and the Meissner effect — both disap- 
pear when the size of a superconducting grain shrinks to a few nanometers jD. Naturally, one 
would like to know if there are remnant effects of bulk superconductivity that differentiate 
normal and superconducting particles in this size range, and if so, how these effects can be 
described at a quantitative level. In a series of experiments on nanometer-sized Al grains, 
Ralph, Black and TinkhamQ have shown that a gap in the excitation spectrum significantly 
larger than the typical level spacing between single-electron eigenstates still exists and is 
destroyed only by a sufficiently strong magnetic field. The size of the energy gap depends 
on the parity of the electron number N on the grain. These results indicate that pairing 
interactions remain important on nanometer scale. 

Several theoretical approaches have been advanced to explain the experimental findings 
mentioned above, and to address the broader issue of superconductivity at nanometer scale 
[3-14]. Attention has been focused on the applicability of the Bardeen-Cooper-Schrieffer 
(BCS) theory|l5| to a situation where (i) the average level spacing becomes comparable to 
the bulk superconducting gap and (ii) the charging energy for adding an electron becomes the 
largest energy scale and effectively fixes N. The converging view is that, under condition (i), 
crossover takes place from the BCS mean-field regime to a fluctuation-dominated regime, 
though remnant pair correlations exist even in the latter case. A general framework for 
addressing the fluctuation effects has been proposed by Matveev and Larking. This is 
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supplemented by numerical calculations and exact solutions H] which offer a quantitative 
account of the crossover. 

Although the problem of small-grain superconductivity under the BCS Hamiltonian is 
essentially solved, one still owes for a simple viewpoint which is generalisable to situations 
where mathematical difficulties prevent an explicit solution. The usual way for achieving 
such a goal is to identify a suitable order parameter which, in the present case, becomes a 
nontrivial task due to condition (ii). Matveev and Larkin proposed to use the ground state 
parity gap Ap = En+i — (En + En+2)/ 1 ^ (N even) of an unpaired electron for this purpose. 
An alternative suggestion was made in Refs.& O, |l2) where a "pair mixing parameter" Amix 
was introduced. Although both quantities carry the effect of pair scattering, they do not probe 
directly the pair coherence which is a defining property of the superconducting state or the 
remnant of it. 

In this paper, we analyse the small-grain behaviour of yet another candidate, the Penrosc- 
Onsager-Yang off-diagonal long-range order (ODLRO) parameter A OD -@, For the BCS 
Hamiltonian, Aod yields directly the strength of a spontaneous symmetry-breaking field, 
whose quantum fluctuations manifest in the ultra-small grain limit. The crossover mentioned 
above is attributed to a gradual weakening of this field against another energy scale of the 
problem, the level spacing S at the Fermi surface. Analytical and numerical results are 
presented to relate Aod to Amix and Ap. 

Let us start with the BCS Hamiltonian for a metallic grain, 

H B cs =J2 ekC kv c k*-9j2 c h c kl Ck 'i Ck 'T- w 

k.cr k,k' 

Here ct (cka) are creation (annihilation) operators of electrons in single-particle eigenstate 
with spin a and energy e&, and g(> 0) is the interaction constant. The sums in Eq. (Q) 
extend over states in the energy range —Tiujd < £fc < Tiojd around the Fermi level, where u>d 
is the Debye frequency. To be definite, we take = kS, where 5 is the average level spacing, 
and k is an integer running from — L to L = hujjj/5. The total number of levels involved, 
SI = 2L + 1, is proportional to the grain volume. The only other dimensionless parameter 
of the model is the relative interaction strength a = g/5, which is expected to depend only 
weakly on grain sizefL^] . 

Under Eq. (Q), levels occupied by a single electron are inert and do not participate in the 
pair scattering process. For the remaining levels, the BCS Hamiltonian can be written in a 
pseudo-spin form]18|| ; 

H BCS =Y / £k(l + 2st)-gY / 4^n (2) 

k k,k' 

where s£ = cj^cj^, = CfejCfc-f, and s| = (c\^Ck] + c^Cfej, — l)/2 are spin-^ operators. In 
terms of the total spin operator S = J^k Sfc ' ^ ne interaction term takes the form, 

H! = -gS+S- = -gS 2 + gS 2 z - gS z . (3) 

Thus the interaction energy £7 is maximised by the state with the largest total spin S but the 
smallest S z (in magnitude), i.e., the state with the largest XY component. On the other hand, 
the single-particle energies = kS, k = —L, . . . . L, act as a nonuniform polarising field in the 
z-direction. Competition of these two tendencies determines the ground state of the system. 

In the usual mean-field analysis, h± = gS^ is treated as a classical spontaneous symmetry- 
breaking field coupled to the XY-component of the pseudo-spins. The strength of this field 
is precisely the BCS gap parameter, i.e., \h±\ = g(S + S* - ) 1 / 2 = Abcs- The phase of S' ± is 
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undefined for a system with a fixed electron number TV which is an eigenstate of 5^[^. As 
the size of the grain shrinks, quantum fluctuations of the amplitude of become important 
when h± no longer dominates over e k in the neighborhood of the Fermi level. 
The Penrose- Onsager- Yang off-diagonal order parameter, 

AoD = ^(H c fcT c H Cfc 'i Cfc 'T>> ( 4 ) 

k,k< 

offers a convenient measure of the strength of the spontaneous symmetry-breaking field in 
both canonical and grand-canonical ensembles. Here (•) denotes either the expectation value 
over a pure state or an ensemble average. It has been shown that, in the thermodynamic limit 
O — * oo, the existence of the ODLRO (in fc-space) is equivalent to Aod being extensive in Q. 
[ [l6[ 0, . The pseudo-spin form of Aod is given by, 

Aod = ^(S + S-) = is XY , (5) 

where Sxy is an effective XY-spin. In the bulk limit, Sxy is also extensive in SI or the grain 
volume. The off-diagonal correlations of the Cooper pairs are reflected in deviations from the 
free-electron results Aod = 1/2 and Sxy = \/^/2. 

We now consider the relation between Aod and the pair- mixing parameter j|, ||, |ll|, p"2] ], 

A mix = J^UfcUfe, (6) 
k 

where 

u k = (cfcxcfcic^c^) 1 / 2 , v k = (4T c fei c fei c fcT> 1/2 : ( 7 ) 

are occupation amplitudes of paired holes and electrons, respectively. In terms of the pseudo- 
spins, we have, 



Amix = £V7-<4> 2 - (8) 

k 

Since pair-scattering yields partially occupied levels, Amix assumes a nonvanishing value under 
([!]). Unlike Aod, however, Amix does not differentiate between coherent and incoherent 
scattering. The following theorem relates Amix to Aod- 

Theorem: For the global ground state of the BCS Hamiltonian ([!]) with P Cooper 
pairs, the following inequalities hold, 

A ^ x <A D<£+A M ixfl-H)- (9) 



Proof: We start with the inequality, 

Aod < ^ + ^ |(4T c fei c fc'l c fc'T>|- ( 10 ) 

An upper bound to each term in the sum is obtained by inserting the identity operator 

I = ^2 h) h\ between the pair creation and annihilation operators and applying the 
Cauchy-Schwartz inequality | J2 n a nK\ 2 < (J2n KI 2 )(E n \ b n\ 2 ), 

|( c feT44 c fe'i c fe'T>| < y ( c h c ki c H c k\) \J (4'T c i'i Cfe 'J. Cfc 'T) 

= v k v k >. (11) 
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Similarly, we can show that 

|(cfc'j.Cfc'T4r c fc4.)| - u k' u k- (12) 
Multiplying (|ll|) with (p|), we obtain, 

K c lT c ii c fe'l c fe'T)| < ^U k V k y/Uk'V kl . (13) 



Substitution of (|l3| ) into (|l0| ) yields 



fc 



The upper bound to Aod in @ is obtained by applying again the Cauchy-Schwartz inequality 
to the sum in (|l4|) with a k = 1 and h k = y/u k v k . 

Unlike the upper bound, the lower bound in (|^) makes specific reference to the BCS 
Hamiltonian for which the interaction energy Ei = — gi7AoD, and to the global ground-state 
"to- Leaving out levels occupied by a single electron, we may write (Q) as, 

u k = <* |l-n fc |*o) 1/2 , v k = <* |«fc|*o) 1/2 , (15) 
where n k is the number operator for a pair. Consider now a normalised BCS-like wave function, 

$ =n(^+^ c i T c ii) io>- (i6) 

The expectation values of -ffecs m the two states and are given by, 

(*o|#bcs|*o) = E K + Ei, (17) 

(*o|#bcs|*o) = E K ~gA 2 Mlx -gJ2<, (18) 

fc 

where Ek — 53k 2efcf^ is the kinetic energy of pairs in the ground state. The lower bound in 
(H) is obtained by noting that ( |l8| ) is an upper bound to (|l|). QED. 

We note in passing that the particle-hole symmetry of -ffp,cs suggests that the global ground 
state (for a fixed parity) is achieved at the half- filling value P — £1/2, though a rigorous proof 
is so far lacking. By the very construction of the BCS Hamiltonian ([!]), the global ground 
state (with respect to all possible values of P) is the one physically realised. 

A direct corollary of the theorem is that, for the ground state of the BCS Hamiltonian, Amix 
is extensive if and only if Aod is extensive. 

In Ref. it was suggested that Amix may become non-extensive when a — g/S is less 
than a certain critical value a c > 0. This is in apparent contradiction with the well-known 
result that the ground state of the BCS Hamiltonian has ODLRO for any a > 0. Since the 
BCS ground state energy -Ebcs provides another upper bound to the true global ground state 
energy Eq, it is possible to show that 

Aod > -Ecpca/dgSl), (19) 
where Ec,bcs is the BCS condensation energy which, for large f2, is given by, 



A 2 

7T, ^BCS 

— J 



1 + 1/1 



, AbCS \2 

- L6 ' 



(20) 



Here Abcs = LS / smh(5 / g) is the BCS gap parameter. Since Ec, bcs is extensive, both Aod 
and Amix are extensive in the bulk limit. 



G.-S. Tian et al.\ Superconducting Correlations. 



5 



Let us now turn to the relations among the spontaneous symmetry breaking field h±, the 
parity parameter Ap, and other spectroscopic parameters. In the bulk limit, \h±\ 3> 8, so 
that pseudo-spins around the Fermi level are completely aligned with the symmetry-breaking 
field. The parity gap Ap is simply the energy cost for removing a level at the Fermi energy. 
This has the effect of reducing Sxy by 1/2 which, from Eq. (||), carries an energy cost 
gSxY = \h±\ — Abcs- Likewise, the energy change for reversing the pseudo-spin, which 
reduces Sxy by one, is 2gSxY = 2Abcs- This is the energy gap for exciting a pair. The 
energy gap for breaking a pair can be derived in a similar fashion and is again given by 2Abcs 
in this limit. It is interesting to explore how these results change for ultra-small grains. 

One possible scheme for addressing pair-scattering effects in the fluctuating regime is a 
renormalised perturbation theory advanced by Matveev and Larking . Treating the interaction 
term as a perturbation and computing the ground state energies to the second order in g, 
Matveev and Larkin obtained, 

A ^f + I ln 2- < 21 > 

The second term in ( [jjj ) can be absorbed in the form of a renormalised coupling constant 
g = g/[l — (g/8) Ln(f2/2)], and finally, Ap can be expressed in a scaling form, 

Ap = A bcs /(<VAbcs), (22) 

where f(x) ~ 1 for x <C 1 and f(x) ~ x/(2lnx) for x ^> 1. Although the scaling ansatz fl2^ ) 
is supported by numerical results, a full-fledged renormalisation group derivation is lacking. 

To check if a similar scaling form applies to the spontaneous symmetry-breaking field h± , we 
have carried out exact diagonalisation of i?BCS for small systems (O = 7, 9, . . . , 21) using the 
Lanczos method. Figure 1(a) shows gSxv / Abcs (diamonds) versus 8/ Abcs for a range of a 
and Q values. For comparison, we have also included data for the parity parameter Ap (open 
circles) and the pair-mixing parameter Amix (solid circles) in a suitably scaled form. The 
scaling ansatz (p2]), which describes well the crossover behaviour of the parity gap Ap, does 
not seem to apply to gSxY and (?Amix in the fluctuation-dominated regime. The lower-bound 
in (|^) (equivalent to Amix < Sxy) ^ s obviously satisfied by the data. Note that, in the small 
grain limit, Sxy is dominated by zero-point fluctuations of rather than the spontaneous 
symmetry-breaking effect. This is reflected in the increase of gSxy from its bulk value Abcs 
seen in Fig. 1(a). 

By computing the ground state wave function to the first order in g, we derived the following 
result, 

AoD-i + f ln2 + 0( ff 2 ). (23) 

Comparison with Eq. ( |2l| ) suggests, on the strongly-fluctuating side, 

A OD * i + (21n2)-p (24) 

Figure 1(b) shows Aod against Ap/8 for different values of a — g/8. At the lower- left corner 
of the graph, there appears to be a good data collapse which compares favorably with Eq. ( |24 ) 
as indicated by the solid line. (Note that part of the descrepancy is due to P/Vt being slightly 
less than 1/2 for odd O.) On the other hand, in the bulk limit, Aqd = ^HAbcs/s) 2 can 
not be expressed as a function of Ap/8 only. This again indicates that the fluctuation and 
correlation effects can not be fully captured by a one-parameter renormalisation group theory. 

In summary, we have established some inequalities that relate the Penrose-Onsager-Yang 
off-diagonal order parameter Aqd to the pair-mixing parameter Amix recently proposed to 
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characterise the crossover from the bulk BCS regime to the fluctuation-dominated regime. It 
is shown that, for the ground state of the BCS Hamiltonian, both quantities are extensive in 
the bulk limit. However, Aod offers a better measure of the superconducting condensation in 
general as it is sensitive also to the "phase" of electron pairs. The nature of the spontaneous 
symmetry breaking for superconducting grains at a fixed electron number N, and relationships 
among various excitations, are examined in the framework of the pseudo-spin representation. 
Our exact numerical results for small systems support the previous view that the fluctuation- 
dominated quantum regime is reached when the BCS gap parameter Abcs becomes comparable 
to the level spacing S. 
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Fig. 1. - (a) The spontaneous symmetry-breaking field gSxy (o), parity parameter (o), and 
pair-mixing parameter Amix (•) against <5/Abcs in suitably scaled form, (b) The ODLRO 
parameter plotted against Ap/S. Solid line indicates result from the perturbative calculation. 
Numbers in legend give the values oia — g/8 for each data set. 



